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Abstract 

The contribution presents a summary of the Gauge/Gravity ap- 
proach to the study of hydrodynamic flow of the quark-gluon plasma 
formed in heavy-ion collisions. Considering the ideal case of a super- 
symmetric Yang-Mills theory for which the AdS / CFT correspondence 
gives a precise form of the Gauge/Gravity duality, the properties of 
the strongly coupled expanding plasma are put in one-to-one corre- 
spondence with the metric of a 5-dimensional black hole moving away 
in the 5th dimension and its deformations consistent with the relevant 
Einstein equations. Several recently studied aspects of this framework 
are recalled and put in perspective. This paper is a written version of 
the four lectures given by the authors on that subject. 
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1 Hydrodynamics are relevant for heavy- ion 
collisions 



One of the most striking lessons one may draw [T, !2] from experiments on 
heavy-ion coUisions at high energy {e.g. at the RHIC accelerator, Brookhaven) 
is that fluid hydrodynamics seems to be relevant for understanding the dy- 
namics of the reaction. Indeed, the elliptic flow [3] describing the anisotropy 
of the low-pT particles produced in a collision at non zero impact parameter 
implies the existence of a collective flow of the particles following a hydro- 
dynamical pressure gradient due to the initial eccentricity in the collision. 
Moreover most hydrodynamical simulations which are successful to describe 
this elliptic flow are consistent with an almost "perfect fluid" behaviour, i.e. 
a small "viscosity over entropy" ratio 7]/s (see, for instance, the reviews [2]). 

The validity of a hydrodynamical description assuming a quasi-perfect 
fluid behaviour has been nicely anticipated in Ref. [4J. The so-called Bjorken 
flow is based on the hypothesis of an intermediate stage of the reaction pro- 
cess, namely a boost-invarian10 quark-gluon plasma phase as a relativistic 
expanding fluid. It is formed after a (quite rapid) thermalization period 
and flnally decays into hadrons, see Fig. [1] The boost-invariance can be 
justifled in the central region of the collision since the observed distribu- 
tion of particles is flat, in agreement with the prediction of hydrodynamical 
boost-invariance, where (space-time) fluid and (energy-momentum) particle 
rapidities are proved to be equal [1], see section 3. 

The Bjorken flow was instrumental for deriving many qualitative and 
even quantitative features of the quark-gluon plasma formation in heavy-ion 
reactions. However, as inherent to the hydrodynamic approach, it says only 
little on the relation with the microscopic gauge fleld theory, i.e. Quantum 
Chromodynamics (QCD). Some important questions remain unsolved, such 
as the reason why the fluid behaves hke a perfect fluid, what is the small 
amount of viscosity it may require, why and how fast thermalization proceeds, 
etc... The problem is made even more difficult by the strong coupling regime 
of QCD which is very probably required, since a perturbative description 

^Thc introduction of hydrodynamics in the description of high-energy hadronic coUi- 
sions has been proposed by Landau 5J, assuming "fuU stopping" initial conditions which 
resuh in a non boost-invariant solution or Landau flow (see 'W for a unified description 
of Bjorken and Landau flows). We will comment later on the relevance of the Landau flow 
for AdS/CFT. 
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Figure 1: Description of QGP formation in heavy ion collisions. The kinematic 
landscape is defined by r = ^y Xq — x'f ; r] = ^ log f^^f^ ; xt = {3^2, 2:3} , where the 
coordinates along the light-cone are XQziz xi, the transverse ones are {x2,X3} and 
r is the proper time, r] the "space-time rapidity" . 



leads in general to a high rj/s. Indeed, the mean free path induced by the 
gauge theory should be small (hence the coupling strong) in order to damp 
the near-by force transversal to the flow, measuring the shear viscosity. 

It is thus interesting to use our modern (but still largely in progress) 
knowledge of non perturbative methods in quantum field theory to fill the 
gap between the macroscopic and microscopic descriptions of the quark-gluon 
plasma produced in heavy-ion collisions. Lattice gauge theory methods are 
very useful to analyze the static properties of the quark-gluon plasma, but 
there are still powerless to describe the plasma in collision. Hence we are led 
to rely upon the new tools offered by the Gauge/Gravity correspondence and 
in particular the one which is the most studied and well-known namely the 
AdS / CFT duality [7] between the Af = 4 supersymmetric Yang-Mills theory 
and the type IIB superstring in the large Nc approximation. The features of 
the gauge theory on the (physical) Minkowski space in 3 + 1 dimensions at 
strong coupling are in one-to-one relation with corresponding ones in the bulk 
of the target space of the 10-d string and in particular in the 5-dimensional 
metric of the AdS space, the boundary of which can be identified with the 
4-dimensional Minkowski space. 

One should be aware when using the AdS/CFT tools that there does 
not yet exist a gravity dual construction for QCD. However, the nice fea- 
ture of the quark-gluon plasma problems is that it is a deconfined phase 
of QCD, characterized by collective degrees of freedom and thus one may 
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expect to get useful information from AdS/CFT duality. This has been al- 
ready proved when describing static geometries by an evaluation of rj/s [S]. 
The subject of the present lectures is the investigation of the Gauge/Gravity 
correspondence, in particular the AdS/CFT duality, in a dynamical setting 
corresponding to a collision. 

2 Relativistic Hydrodynamics and Bjorken Flow 

On theoretical grounds, there are quite appealing features for applying hy- 
drodynamic concepts to high-energy heavy-ion reactions. Such concepts have 
been already introduced some time ago 0, S] and find a plausible realization 
nowadays. The fact that a rather dense interacting medium is created in the 
first stage of the collision allows one to admit that the individual partonic or 
hadronic degrees of freedom are not relevant during the early evolution of the 
medium and justifies its treatment as a fluid. For the same reason local equi- 
librium is a plausible assumption. Moreover, the high quantum occupation 
numbers allow one to use a classical picture and to assume that the "pieces 
of fluid" may follow quasi-classical trajectories in space-time, expressed as 
an in-out cascade [9] with straight-line trajectories starting at the origin (see 
Fig. 2), with 

y = v (1) 

where 

are respectively the rapidity and "space-time rapidity" of the piece of the 
fiuici. 



Note that ([T]) can be rewritten in the form 

2y = logM^ — logu~ = logx^ — logx~ (3) 

where = e^^ are the light-cone components of the fluid (four-)velocity 
and = xo± xi are the light-cone kinematical variables. 

^We keep the conventional notation ?/, not to be confused with viscosity. The difference 
is clear enough to avoid mistakes. 
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Figure 2: In- Out cascade. The "piece of fluid" with space-time rapidity 
rj gives rise to hadrons at rapidity y = t], after crossing the "freeze-out" 
hyperbola at fixed proper-time r. 

Taking ([T]) as the starting point and using the perfect fluid hydrodynam- 
ics, Bjorken developped in his seminal paper |4j a suggestive (and very useful 
in many applications) physical picture of the central rapidity region of highly 
relativistic collisions of heavy ions. In this picture the condition ([1]) leads to 
a boost-invariant geometry of the expanding fluid and thus to the central 
plateau in the distribution of particles. 

Let us introduce the relativistic hydrodynamic equations in light-cone 
variables. We consider the "perfect fluid" approximation for which the 
energy-momentum tensor is 

T"'^ = (e + - pri"" (4) 

where e is the energy density, p is the pressure and u'^ is the 4-velocity. We 
assume that the energy density and pressure are related by the equation of 
state: 

e = gp (5) 



5 



where l/y^ is the sound velocity in the hquid. For the "conformal case" 
T^''' = and thus ^ = 3. 
Using 

= M° ± = (6) 

and introducing 

,._o±,._,.,^(^±^,^,^, (7) 



where r = Vx+x is the proper time and t] is the spatial rapidity of the 
fluid, the hydrodynamic equations 

d^.Tf'" = (8) 

take the form 

5^T°^ + ^d+{T'^ ± T°°) - ^9_(r" T = . (9) 
Using now (jlj) and the equation of state ([5]) we deduce from this 
9d+[\ogp] = — d+y —e ^d^y 

gd_[\ogp] = ^l±^d^y + ^^e^yd+y . (10) 

These are two equations for two unknowns which describe the state of the 
liquid: the pressure p and the rapidity y. They should be expressed in terms 
of the positions a;"*", x~ in the liquid. Other thermodynamic quantities can be 
obtained from the equation of state ([5]) and the standard thermodynamical 
identities: 

p + e = Ts] de = Tds (11) 

where we have assumed for simplicity vanishing chemical potential. 
The result is 

e = gp = eoT3+^ ; s = SqT^ s ~ e^/^^+i). (12) 

The simplest possibility to describe the expansion of the fluid was sug- 
gested by Bjorken [4j who proposed to use the Ansatz ([T]) in the hydrody- 
namical context. Introducing ([T]) into ( JTOl) we obtain 

gd+[\ogp] = -^^^ ; gd-[\ogp] = -^^J^ (13) 
6 



from which we deduce 

where po is a constant, and thus specifically 



p = e/3 = po (a:+x-)-2/3 = r-^/=^ oc (15) 

for the conformal case. 

Thus the system is boost-invariant: the pressure does not depend neither 
on T) nor on y. So are e, s and T, given by (IT^ . It is interesting to note that 
the Landau flow corresponds asymptotically only to a logarithmic correction 
of relation ([1]), namely 

~ x^^/log x"^ , (16) 

which gives flnally rise (as already noticed in [5], and for instance recently 
discussed in [10]) to a gaussian shape in the y distribution of the entropy, 
revealing a non boost-invariant picture, at least at some distance from central 
rapidity. 



3 Interest of AdS/QCD duality 

In the previous sections we mentionned the ubiquity of hydrodynamic meth- 
ods in the description of QGP produced at RHIC. Yet, despite their success 
in describing data, we have to keep in mind that they are used as a phe- 
nomenological model without a real derivation from gauge theory. This is 
quite understandable since almost perfect fluid hydrodynamics is intrinsi- 
cally a strong coupling phenomenon - for which one lacks a purely gauge 
theoretical methoco. 

On the other hand there exists a wide class of gauge theories, which 
can be studied analytically at strong coupling. These are super conformal 
field theories with gravity duals. String theory methods (namely AdS/CFT 
correspondence) maps gauge theory dynamics (CFT) at strong coupling and 
large number of colors into solving Einstein equations in asymptotically anti- 
de Sitter space (AdS). The theories with gravity duals can differ substantially 

•^Lattice QCD methods do not work well here as this would require analytical continu- 
ation to Minkowski signature which is nontrivial in this context 
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from real world QCD at zero temperature. The best known example of 
such theory - A/" = 4 super Yang-Mills (SYM) - is a superconformal field 
theory with matter in the adjoint representation of the gauge group SU(Nc). 
Because of the conformal symmetry at the quantum level this theory does 
not exhibit confinement. On the other hand differences between Af = 4 SYM 
and QCD are less significant above QCD's critical temperature, when quarks 
and gluons are in the deconfined phase. Moreover it was observed on the 
lattice that QCD exhibits a quasi-conformal window in the certain range of 
temperatures, where the equation of state is well-approximated by e = 3p. 
The above observations together with experimental results suggesting that 
quark-gluon plasma is a strongly coupled medium is an incentive to use the 
AdS/CFT correspondence as a tool to get insight into the non-perturbative 
dynamics. 

4 AdS/CFT setup 

We will now describe how to set up an AdS/CFT computation for deter- 
mining the spacetime behaviour of the energy- momentum tensor [11]. This 
method does not make any underlying assumptions about local equilibrium 
or hydrodynamical behavior. We will obtain hydrodynamic expansion as a 
generic late time behaviour of the expanding strongly coupled plasma. 

Suppose that we consider some macroscopic state of the plasma charac- 
terized by a spacetime profile of the energy-momentum tensor 

T,.{xn . (17) 

Then, since the AdS/CFT correspondence asserts the exact equivalence of 
gauge and string theory, such a state should have its counterpart on the string 
side of the correspondence. Typically it will be given by a modification of 
the geometry of the original AdS^ x metric. This comes from the fact 
that operators in gauge theory correspond to fields in supergravity (or string 
theory). When we consider a state with a nonzero expectation value of 
an operator, the dual gravity background will have the corresponding field 
modified from its 'vacuum' AdS^ x value. In the case of the energy 
momentum tensor the corresponding field is just the 5-dimensional metric. 
One then has to assume that the geometry is well defined i.e. it does not 
have a naked singularity - a singularity not hidden by an event horizon - 
. This principle will select the allowed physical spacetime profiles of gauge 
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theory energy-momentum tensor. Thus together with the Einstein equations 
this becomes the main dynamical mechanism for the strongly coupled gauge 
theory. 

The simplest way to formulate the precise correspondence between the 
expectation value of the energy-momentum tensor and bulk geometry is to 
use the Fefferman-Graham system of coordinates [12] for the latter: 



9iJ.u{^'^ 1 z)dx^dx + dz (18) 

This metric has to be a solution of 5-dimensional Einstein's equation with 
negative cosmological constant^: 

R^lu - ^9^luR -Q9f,u = . (19) 

The expectation value of the energy momentum tensor may be easily re- 
covered by expanding the metric near the boundary z = 0, following the 
"holographic renormalization" procedure [13], 



g^^ix'', z) = + z'g^^J (x") + . . . . (20) 

Then 

{T,.{xn) = ■ gj^K^n ■ (21) 

This relation can be used in two ways. Firstly, given a solution of Ein- 
stein equations we may read off the corresponding gauge theoretical energy- 
momentum tensor. Secondly, given a traceless and conserved energy-momentum 
profile one may integrate Einstein equations into the bulk in order to obtain 
the dual geometr}!^!. Then the criterion of nonsingularity of the geometry 
obtained in this way will determine the allowed spacetime evolution of the 
plasma. Let us note that this formulation is in fact quite far away from a 
conventional initial value problem. 

Before we move to the case of expanding plasma, it is convenient to con- 
sider the simple situation of a static uniform plasma with a constant energy 



^One can show that such solutions Uft to lO-dimensional solutions of ten dimensional 
type IIB supergravity. The effective 5-dimensional negative cosmological constant comes 
from the 5-form field in 10-dimensional supergravity. 

^This can be done order by order in z^, which is a near-boundary expansion. However 
potential singularities are hidden deep in the bulk, thus this power series needs to be 
resummed. 



9 



momentum tensor. Then the Einstein's equations can be solved analytically 
and we find [TT] that the exact dual geometry of such a system is 



"^^ - - (TTW *^ ^ + '''"^^ ^ ^ ■ '''' 



This metric may look at first glance unfamiliar, but a change of coordinates 

z 



1 + 

^0 



(23) 



transforms it to the standard AdS Schwarzschild static black hole 



ds' = -^-^P^dt- + ^ + -4^^ (24) 

Z-^ Z^ 1 — Z^/Zq 

with zq = zq/\/2 being the location of the horizon. Before we proceed further, 
let us note here one crucial thing: the fact, that the dual geometry of a 
gauge theory system with constant energy density is a black hole was not an 
assumption, but rather an outcome of a computation. 
The Hawking temperature 

T = ^^^ (25) 

nzo Tizo 

is then identified with the gauge theory temperature, and the entropy with 
the Bekenstein-Hawking black hole entropy 



Nl_ TT^ 



which is 3/4 of the entropy at zero coupling. To finish our discussion of the 
static black hole, we note that the Fefferman-Graham coordinates cover only 
the part of spacetime lying outside the horizon. 



5 Boost invariant flow 

Let us now apply the above procedure to a generic boost-invariant fiow, in 
view of making contact with the hydrodynamical Bjorken fiow described in 
section 2. However we do not want to make any preassumptions on the 
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dynamics, since we would like to recover the hydrodynamic behaviour as an 
outcome of an AdS/CFT computation. To this end let us consider the most 
general gauge theory energy-momentum tensor which is boost-invariant and 
does not depend on transverse coordinates (see Fig. 1). Then conservation 
of energy-momentum d^T^'^ = and tracelessness = allow to express 
all nonvanishing components of T^^, in terms of a single function £:(r) - the 
energy density at mid rapidity: 



T : 



( sir) \ 

-r^j-e{r)-T\{T) 

e(r) + ir|:e(r) 

e[T)^\Tj-^e[T) I 



(27) 

Let us concentrate, following [11] on the late time asymptotics of this 
function i.e. 

e(r) ~ ^ + . . . (28) 

for r ^ oo. Energy positivity requires that < s < 4. We will consider sharp 
inequalities herqj. The most general metric consistent with the symmetry 
assumptions is 

ds^ = ^ + — . 29 

z^ 

In order to find the late time form of the solution corresponding to sir) = 
1/r* we may solve the Einstein equations in a power series for the metric 
coefficients 



N 



[r)2^+'" (30) 

n=0 

where ao(''") = ~£^(t) = — l/r''. Then from each coefficient a,i(r) we may 
extract the leading large r behaviour and neglect the subleading terms. It 
turns out that this procedure is exactly equivalent to introducing a scaling 
variable 

v = ^ (31) 

and assuming the metric coefficients to be just functions of v e.g. a{z, r) = 
a{v) in the large proper time limit (namely r — >■ oo, z ^ oo with v kept fixed). 



^Recently the case s = 4 has been considered in [T4] 
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In this limit Einstein's equations become just ordinary differential equations 
and may be solved analytically. The singularity of these geometries can then 
be tested by computing the scalar curvature invariant 

^2 ^ Rf^--PR^^^^ . (32) 

Since our solutions are defined only in the large proper time limit r oo 
with the scaling variable v kept fixed, we have to evaluate 9^^ in the same 
mannei0. 

This procedure is described in detail in [TT]. The result is that 

• for generic s the resulting solution is singular 

• the only nonsingular solution corresponds to s = | which is just the 
hydrodynamic Bjorken expansion (see f|T5l) . section 2) 

• the resulting metric takes the form 



1 _ eg z 

^ 3 r4/3 



1 + £0. 



1473 



1 _l_ ep z 

3 r4/3 



(r^rfy^ + dx]_ 



+ 



dz' 



(33) 



where we reinstated the dimensionful parameter cq so that 

£(r) = eo/rt . (34) 

Let us note some salient features of this result. The geometry (!33l) bears 
striking resemblance to the AdS black hole geometry (!22l) but with the posi- 
tion of the 'effective horizon' being time dependent 




(35) 



^It should be stressed however, that this condition is reaUy a condition of regularity 
of the expansion of the curvature invariant. It is safe to do as long as each term in the 
large proper-time expansion is regular. On the other hand any singularity present in this 
expansion might be either a genuine curvature singularity or a singularity of the expansion, 
see a detailed discussion in section [H 
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Then assuming similar relations as for the black hole case one gets the Bjorken 
scaling of the temperature and entropy. 

T = ^ = Hr-^ 

S OC ^ = const . (36) 

We see that the 'movement' of the horizon into the bulk of AdS corresponds 
physically to cooling of the expanding gauge theory plasma system. 

A significant fact that has to be kept in mind is that the geometry (133|) . 
in contrast to (!22l) . is not an exact solution of Einstein's equation. It is valid 
only for large times. For smaller times it has to be modified. We will now 
discuss this issue in more detail as it reflects important physical properties 
of the gauge theory plasma. 



6 Beyond perfect fluid 



The geometry fl33|) is only a solution of Einstein equations in the scaling limit. 
However with some effort, one can get also the first sub leading corrections 
to the metric i.e. ^ 

a{z, t) = ao{v) H — j a2{v) + . . . (37) 

T3 

Then after evaluating IR^, keeping track of sub leading terms we find 

= Ro{v) + ^R2{v) + ... (38) 

T3 

where -Ro(^^) is finite, but R2{v) develops a 4*^ order pole singularity. The 
physical meaning of this behaviour is indeed quite clear. The geometry (1371) 
is dual to a state in gauge theory which undergoes expansion according to 
exact perfect fluid hydrodynamics. Yet we know that gauge theory plasma 
has nonzero viscosity and hence the perfect fluid behaviour 

e{r) = 4 (39) 

T3 

is not exact but, if it would be described by viscous Bjorken expansion (vis- 
cous hydrodynamics), it would be modified to 

.(r) = 4fl-^ + ...) (40) 

T3 \ T3 I 
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where tiq is related to the shear viscosity through t] = tiq/t (which follows 
from the scaling 77 oc T^). 

Let us show how this arises using the AdS/CFT methods. We will not 
presuppose a specific form of subleading correction but will start from 

.M^^(i-f?....) («) 

with a generic r. In order to verify that plasma expansion follows viscous 
hydrodynamics we will have to first show that r = ^. The metric coefficients 
will now have an additional piece scaling ). It turns out that the 

curvature scalar 9^^ is always nonsingular at that ordei@. Hence we have 
to go one order further i.e. find all coefficients appearing in the following 
expansion 

a{z, t) = ao{v) + ^ar{v) + -^a2r{v) + ■\a2{v) + ... (42) 



Then the curvature scalar has the form 

T T 7-3 



= R^(y) + -R,{v) + —^R2r{v) + —R2{v) + ... (43) 



with Roiv) and Rriv) being nonsingular, while both R2r{v) and R2{v) turn 
out to have 4*^^ order pole singularities. In order for them to have a chance 
to cancel we have to have 

r = ^ (44) 

which is exactly the scaling of a viscosity correction to Bjorken flow. More- 
over cancelation occurs only when the shear viscosity coefficient has the 
valu^ 

r/o = 2-^3-t (45) 

which is equivalent to rj/s = I/Atc (for details see [IS]). In a similar man- 
ner one can go one order higher and determine a coefficient of second order 
hydrodynamics. However at that order, it turns out that there remains a 
leftover logarithmic singularity. We will show, in section [HI that the logarith- 
mic singularity arises due to a pathology of the Feffer man- Graham expansion 
and can be avoided when one makes a different late time expansion. 



^This was first observed for r = 2/3 in [15] . 
^We set liere eo = 1. 
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Finally let us comment on why it is interesting to verify the viscous hy- 
drodynamic behaviour with the specific viscosity coefficient for the expand- 
ing plasma. Already before, there have been studies of linearized perturba- 
tions around the uniform plasma which demonstrated that hydrodynamic 
behaviour appears for small ffuctuations and the value of viscosity was ob- 
tained from the Kubo formula [S]. It was interesting to verify whether hy- 
drodynamics also applies in its fully nonlinear regime. The agreement of 
the resulting value of the viscosity coefficient is thus a nontrivial consistency 
check. 

Another motivation for developing an AdS / CFT framework for studying 
such time-dependent phenomenae is the fact that some of the most interesting 
and puzzling phenomena in heavy ion collisions are definitely very far from 
equilibruum. We will mention some examples in section [91 



7 Beyond boost-invariance 

The calculations presented in the previous sections were performed for sys- 
tems with boost invariance symmetry and full translational and rotational 
symmetry in the transverse plane. This allowed us to perform explicit com- 
putations as the symmetry assumptions effectively reduced the calculation 
to systems of ordinary differential equations. In this manner we obtained 
directly the solution for gauge theory energy density ^(t). Then, in order to 
find the link with hydrodynamics, we found that this solution is a solution 
of hydrodynamic equations with specific values for the transport coefficients. 

This approach has both an advantage and a drawback. The advantage is 
that one does not presuppose any kind of dynamics and one may try to apply 
it in contexts very far from equillibrum, where hydrodynamic description does 
not apply. The drawback is that the appearance of hydrodynamic equations 
is not transparent and it is difficult to relax the symmetry assumptions due 
to the complexity of solving nonlinear Einstein's equations. 

Recently the latter drawback was addressed and it was shown in general 
how the equations of hydrodynamics arise from the gravity side [17]. Here 
we will briefiy review this approach. 

Let us start from the static black hole fl24ll22p but written in yet another 
coordinate system - the incoming Eddington-Finkelstein coordinates: 

ds"^ = -2dtdr -r^ (l- df + r\jdx'dx^ . (46) 
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Here T is the temperature, r = oo corresponds to the boundary, x'^ = const 
are null curves going from the boundary into the black hole. The advantage 
of this coordinate system is that it is well defined on the horizon and extends 
all the way from the boundary to the singularity at the center of the black 
hole. 

The geometry given above corresponds to a uniform plasma at rest (i.e. 
with the 4-velocity = (1, 0, 0, 0)) and given temperature T. We may now 
perform a boost (and perform a dilatation) to obtain the dual geometry to 
a moving plasma system with uniform 4-velocity u'^ and temperature T: 



ds"^ = —2u^dx^dr — r'^ I 1 ^ | u^Uydx^dx^ + r^{r]^j, + u^Ujj)dx^dx^ (47) 



The idea of ref.[T7] is to allow and T to be (slowly- varying) functions of the 
spacetime coordinates. Once this is done the geometry fH7|) ceases to be an 
exact solution of Einstein equation because of nonvanishing gradients of the 
parameters and T. This suggests to perform an expansion of the solution 
in terms of gradients which has been carried out in [17] up to second order in 
derivatives. The integration constants arising at each order are again fixed 
by requiring regularity of the metric at the horizon. The resulting metric is 
expressed in terms of 4-velocities and temperatures and their derivatives, so 
when one extracts the energy-momentum tensor it will be given directly in 
terms of those quantities. Up to first order the expression is 



The first term is just the perfect fluid energy momentum tensor, while the 
second term involves the shear viscosity. This result essentially demonstrates 
how general hydrodynamic equations arise from gravity in AdS/CFT. Indeed, 
once it is shown that the general form of the gauge theory energy-momentum 
tensor has the form ( HHl) . then conservation of energy momentum d^T^^ = 
is equivalent, by definition, to conformal relativistic Navier-Stokes equations. 
As a byproduct, the above construction also gives a map from solutions of 
(viscous) hydrodynamics to gravity solutions. However this setup requires 
that the starting point is not far off from equilibruum. For processes which 
do not admit a hydrodynamic description (like the early stage of a heavy-ion 
collision) one has to resort to different methods. 





shear J 



(48) 
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8 Reduction of Singularities 



The leftover logarithmic singularity found in the third order of the square of 
the Riemann tensor [18] (as well as in the higher curvature invariants [T9] ) 
might be the signal of either genuine curvature singularity or the singularity of 
the chosen expansion schemq^j. If the first is true, this means that the whole 
framework is inconsistent and either one needs to include additional degrees 
of freedom to cure it or the boost-invariant flow is unphysical"1. The results 
presented in [19] show that no supergravity field can fix the problem, which 
led to conjectures, that boost-invariant flow cannot be realized within the 
supergravity framework [20] . On the other hand, the gravity dual of general 
fluid flow up to the second order in derivatives was shown to be regular and 
it was hard to imagine how possible singularities could form in the third 
order [TTl [21]. The resolution of this puzzle was presented in [22] (see also 
|23j). It turns out that there exists a singular coordinate transformation from 
Fefferman-Graham coordinates to Eddington-Finkelstein ones, which yields a 
completely regular and smooth metric from the boundary up to the standard 
black-brane singularity. This leads to regular (apart from the standard black- 
brane singularity) large proper-time expansion of curvature invariants. The 
metric ansatz in Eddington-Finelstein coordinates takes the form 



ds^ = 2df ■ dr - r^A (f , r) df ' + (1 + r ■ ff e^^^'^'My^ + r^e^'-^'^'Mxl (49) 

and was motivated by the boosted black-brane metric fH7|) with a boost and 
dilatation parameters u = 1-df and T ~ f~^/^. The functions A (f, r), b (f, r) 
and c (f, r) are expanded in a large-proper time expansion analogously as it 
was in the Fefferman Graham case, i. e. 



A (f , r) = io (r • f'/') + (r ■ f'/'') + ^A^ (r ■ f'/') +... (50) 

This form of expansion can also be justifled by [23]. The terms damped by 
inverse power of proper time come from the gradient expansion. The bound- 
ary metric in proper-time-rapidity coordinates has non-vanishing Christoffel 

^*^As it was stressed before, the large-proper time expansion of curvature invariants is 
not diffeomorphism-invariant. The encountered singularities are physical only if there is 
no coordinate transformation which removes them. 

Since it corresponds to the naked singularity on the gravity side 
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symbols F ~ f~^, thus the four velocity gradient Vm (which is constant in 
these coordinates) gives a factor of f~^. On the other hand the expansion 
parameter multiplying each term in gradient expansion is the inverse power 
of the temperature T(see Because T ~ f~^^^, the overall damping is 

indeed fact derived in [18] from the non- singularity argument. 

The non-perturbative0 piece in the metric at dy^ introduced in [231 
responsible for a correct limit energy density -^0. It also becomes impor- 
tant if one wants to solve the problem of early-time dynamics [27] using 
Eddington-Finkelstein coordinates. 

The integration constant^ are fixed by requiring the regularity of the 
metric functions Ai (v), hi {v) and Cj {v) at each order i. This is justified since 
the Eddington-Finkelstein are valid for f > and Q < v = r ■ f < oo. The 
singular coordinate transformation which takes the metric from Eddington- 
Finkelstein coordinates to Fefferman-Graham ones is given order by order in 
the gradient expansion by 



r (r, z) 
r (r, z) 



T ■ 

1 

Z 



Z ■ T 



-1/3 



) + Z^T, ( 



Z ■ T 



-1/3 



-V3) +...}, (51) 

) + ;^i?,(z-r-/3)+...}. (52) 



The leading-order solutions (corresponding to the perfect fluid on the gauge 
theory side) are related by 



r 



r 1 



^2/3 



31^ 31^ 

31/4 



tan 



31^ 
V2 



r ■ T 



1/3 



+ 



+ 



4^/2 



log 



r • r 



1/3 V2 ' 



1 

z 



3-r4/3' 



(53) 



The transformation is singular at z = 3^^^t^^^, which is precisely the locus 
of the logarithmic singularity encountered in [18] . Formulas for higher order 



^^In the sense of large-proper time expansion 

^•^Not all of them - there is a remaining gauge freedom (coordinate transformation), 
which leaves the metric ansatz unchanged: r r + f{f), where /(f) is an arbitrary 
function 
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transformation coefficients are too long to be presented here and can be 
found in The energy-momentum tensor extracted from the solution in 
Eddington-Finkelstein coordinates reproduces the energy momentum tensor 
obtained in |18]. 

9 Beyond hydrodynamics 

Gauge-gravity duality has already proven to be an invaluable tool in de- 
scribing properties of static or near-equilibrium (hydrodynamics) strongly 
coupled gauge theory systems. Noticeable achievements in that direction are 
the viscosity evaluation bound [S] and the consistent formulation of the sec- 
ond order conformal hydrodynamics ^221 [IZj. These successes came from the 
holographic understanding of hydrodynamics. On the other hand there is 
much more interesting and nontrivial dynamics than hydro. Far from equi- 
librium behavior of gauge theories is a fascinating and pretty much open 
problem of experimental importance, like the early universe or initial stages 
of heavy ion collision^. The AdS/CFT correspondence is surely capable to 
shed new light on these problems, or even be understood as a formulation of 
far from equilibrium gauge theory. 

In the context of heavy-ion collisions the most important and probably 
the most difficult questions concern the issues of early time dynamics ^27] 
and the transition to an isotropic [31] and thermalized medium. One of 
the puzzles here is the short time in which nuclear matter approach local 
equilibrium. Experimental data fitting well with hydrodynamical simulations 
with small viscosity justified applications of the AdS/CFT correspondence 
at strong coupling for the late stages of heavy-ion collisions. It is not clear 
to what extent early time dynamics is driven by non-perturbative effects and 
whether the lessons learned from AdS / CFT might be directly applied to the 
nuclear matter in the early stages of the evolution. Approaching equilibrium 
is also of an interest from the General Relativity point of view. Isotropic 
and thermalized matter on the gauge theory side corresponds to a black hole 
in AdS, whereas thermalization and approach to local equilibrium should be 

^"^In the late stages of heavy ion colhsion, strongly coupled quark-gluon plasma forms 
and holographic technics at strong gauge coupling are better justified then just after the 
collision (running of the coupling) . Nevertheless applying AdS / CFT correspondence to 
describe far from equilibrium processes in gauge theories is an interesting problem even 
from a purely theoretical point of view. 
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governed by the dynamics of gravitational collapse. 

Perhaps some of these questions might be answered by studying collisions 
of shock-waves in AdS. The geometry corresponding to a projectile in 3+1 
dimensions was constructed in [11] using holographic renormalization. The 
metric 

ds^ = ^ |-2da:+dx- + / (a;") ■ (dx")^ + dxl + dz^j (54) 

with an arbitrary function / (x~) corresponds to the situation when 

• the dynamics is one- dimensional (i.e. no dependence on transverse 
coordinates), 

• the energy-momentum tensor depends only on a single light-cone vari- 
able (here chosen to be x~ = x° — x^). 

Traceless and conserved energy momentum tensor satisfying the above 
assumptions takes a particularly simple form - its only non-zero component is 
T = / (x^). Choosing / (x^) = M6 (x^) leads to a shock- wave - infinitely 
thin plane of matter moving at the speed of light, which is a toy- model 
for highly boosted nucleus. The idea is to collide two such projectiles and 
single out the physical behavior of the plasma from the regularity of the dual 
geometry. This is a difficult problem, because of the broken boost-invariance, 
which leads to solving Einstein equations in 3 variables (x"*", x~ and z or 
equivalently r, y and z). The geometry before the collision (x"*" + x~ < 0) is 
known - it is simply the superposition of two incoming shock-waves 

ds^ = \ |-2dx+dx- + M5 (x-) z^ {dx-y + MS (x+) z^ {dx-^f + dxl + dz^^ . 

(55) 

Shock-waves collide at x"*" = x =0 and from now on the dynamics of 
the system must be deduced from Einstein equations. The first attempt to 
address this issue in [31] focused on the simpler setup than presented so far - 
a shock- waves collision in l-fl dimensions. The energy- momentum tensor for 

such a system before the collision is given by T++ = / (x~) and T = g (x^) 

with vanishing off-diagonal components. This is at the same time the most 
general form of the energy- momentum tensor for a 1 + 1 dimensional CFT. 
A nice feature is that the dual geometry for the whole collision process can 



20 



be constructed here exactly. However in this low dimensional context, the 
projectiles pass each other unaffected (or propelled back-to-back [21]), so 
the physics of plasma production and thermalization is absent here. On 
the other hand the problem of genuine interest - collision of shock-waves 
in 3+1 dimensions - requires some approximation scheme in which Einstein 
equations become tractable. Up to now, two proposals have been made. 
The first one [52] treats proper time as a small parameter but suffers from a 
negative energy density in some regions due to the conditions imposed at the 
light-cone. The second one [33] solves Einstein equations perturbatively in 
M leading to the prediction that shock-waves stop almost immediately after 
the collision (reminding of the full stopping condition of the Landau flow, see 
section [21). This problems surely deserve further studies. 

There are also other studies of dynamical processes in an evolving plasma 
system which go beyond hydrodynamics. One example is the problem of 
thermalization of small perturbations around the expanding plasma (some 
first investigations has been performed in [H]). Another use of the evolving 
geometries is to study other physical processes in the presence of the evolving 
plasma system like e.g. the physics of mesons and flavours studied through 
embedding D7 branes in the time-dependent geometries [35j. Finally one may 
study isotropisation of anisotropic plasma. The first investigations have been 
performed in [3l] (see the talk by P. Witaszczyk in the same proceedings). 

10 Summary 

The Gauge/ Gravity approach to the formation and evolution of a quark- 
gluon plasma in heavy-ion collisions described above has the interest of cast- 
ing an exploratory bridge between the rigorous results of string theory and 
some pending questions raised by the experiments on quark-gluon plasma. 
These questions cannot yet be raised in the framework of strongly coupled 
QCD, for which we do not possess the adequate tools, but they can be ad- 
dressed for the first time in a quantitative and rigorous way in the supersym- 
metric case of the AdS/CFT correspondence. It is thus a novel and valuable 
approach and can serve as a model for further studies. 

Let us summarize some aspects of this approach, being aware (and with 
apologies for those not quoted or mentioned) that this subject is in constant 
development which will force us to mention only a few of them. 

Starting with the experimental evidence that hydrodynamics is relevant 
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in the formation and evolution of a quark-gluon plasma in heavy-ion col- 
lisions and in particular of the "Bjorken flow" description, we present the 
AdS/CFT setup allowing to describe the dynamics of the plasma (in the 
AdS/CFT case). We show that it is possible to derive the geometry dual to 
the asymptotic evolution of the plasma in terms of an expansion in a scal- 
ing variable, the nonsingularity requirement on the Gravity side gives a set 
of "selection rules" on the Gauge side: the perfect fluid at first order, the 
T}/ s = I/Att property and other transport coefficients at higher orders... 

Beyond the boost-invariant Bjorken flow, there exists an intriguing but 
rigorous one-to-one correspondence between the complete (and even com- 
pleted using AdS/CFT!) hydrodynamic equations and the solutions of the 
Einstein equations in the bulk of the 5-dimensional space. Some apparent 
obstacles, such as the appearance of logarithmic singularities in the asymp- 
totic expansion of the 5-dimensional metric, have been shown to be a mere 
artefact of the choice of the expansion parameter and have been cured. 

Beyond the hydrodynamical description of the transient plasma, one goal 
is now to exoplore the dynamical aspects far from equilibrium. We describe 
some very recent attempts, which even though not conclusive yet, show the 
interest of the extension of Gauge/Gravity correspondence to attack some 
down-to-earth problems, such as the short thermalization time, probably 
observed by the heavy-ion phenomenology and more generally the effect of 
the initial conditions on the whole process. 

It is quite interesting to see that some complex aspects of Gauge field 
theory dynamics can find unexpected answers from Gravity. It would be 
intringuing that some nontrivial aspects of Gravity (such as the dynamics of 
moving black holes) also could gain some new insight from the correspondence 
with some aspects of heavy-ion collisions. 
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